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Interfacial friction plays a crucial role in the mechanical properties of carbon nanotube based 
fibers, composites, and devices. Here we use molecular dynamics simulation to investigate the 
pressure effect on the friction within carbon nanotube bundles. It reveals that the intertube frictional 
force can be increased by a factor of 1.5^4, depending on tube chirality and radius, when all tubes 
collapse above a critical pressure and when the bundle remains collapsed with unloading down to 
atmospheric pressure. Furthermore, the overall cross-sectional area also decreases significantly for 
the collapsed structure, making the bundle stronger. Our study suggests a new and efficient way to 
reinforce nanotube fibers, possibly stronger than carbon fibers, for usage at ambient conditions. 
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Carbon nanotubes (CNTs) are considered the 
strongest and ideal reinforcing fibers due to their ex- 
ceptional mechanical properties, low density, and high 
aspect ratio. However, although the axial strength and 
stiffness of individual CNTs are of the order of 50^100 
GPa and 1 TPa, respectively [1-5], the highest strength 
of CNT and CNT-reinforced fibers, ranging from 0.85 to 
3.3 GPa [Mil, is nearly 2^3 orders of magnitude lower 
than individual tubes and about 1/3^1/2 of the strongest 
Toray carbon fibers [l3[ • As good alignment improves the 
translation of axial properties of individual tubes to those 
of the fiber, efforts were reported to grow ultr along and 
well- aligned CNT arrays (forests) |7[, and to improve the 
direct spinning method [10[. Post-spin treatments, e.g., 
infiltration, twisting, heating, and stretching, have been 
reported to improve the load transfer between CNT bun- 
dles, by making better oriented network [ll| and closed 
packing [12[ of bundles. However, CNTs often do not 
exist as individual tubes but group into bundles, the ba- 
sic component of the spun fibers. Therefore it can be 
of great importance to improve the strength of bundles. 
So far, the bundle strength is reported to be about 10 
GPa for the length of several micrometers [3]. Consider- 
ing the strength loss from component filaments to tradi- 
tional fibers [141 ]. it is hard to achieve the same strength 
as carbon fibers by grouping CNT bundles through var- 
ious spinning treatments. One problem strongly related 
to the bundle strength is that individual tubes in bun- 
dle tend to slide easily against each other [15[ . Recently, 
translational static and sliding frictions in a multi-wall 
CNT were measured to be 0.014 and 0.009 me v/A per 
atom, respectively [16| , and the frictional force in a bun- 
dle was reported, surprisingly, of several orders of magni- 
tude greater [l?], EH , due to different experimental con- 
ditions and probably also the existence of impurities. Al- 
though still hard to measure the friction between defect- 
free tubes, there is no doubt that in order to achieve a 
strong bundle, tubes within it should be sufficiently long. 

Here we show, from a series of molecular dynamics sim- 



ulations, that pre-pressing on CNT bundles can greatly 
enhance intertube friction and consequently the strength 
of the nanotube fibers. The underlying mechanism in- 
volves the structural transition of CNTs, accompanied 
by an increase of intertube frictional force and the de- 
crease of the cross-sectional area as well. All tubes col- 
lapse above a critical pressure and remain collapsed af- 
ter unloading, especially the large tubes. The friction 
increase, by a factor ranging from 1.5 to 4, is chirality 
dependent, and is strongest for nonchiral tubes. Taking 
into account this new feature of pre-pressing, it might 
be possible to spin nanotube fibers stronger than carbon 
fibers under current spinning techniques. Furthermore, 
although strongly related to previous studies on radial 
mechanical translation during the transition of individ- 
ual CNTs or CNT bundles \±<M2^, our study focusing on 
axial translation is obviously new and shows inspirational 
results. 



RESULTS AND DISCUSSION 

The simulation is set for defect-free single-wall CNT 
bundles which, named by tube number in the box 
and tube chirality, are 16 x (23,0), 16 x (40,0), and 
16xrandom-chirality (RC), respectively. [T] shows the 
volume-pressure (V-P) relation during the pressure load- 
ing up to 1.2 GPa and unloading down to GPa for 
the 16x(23,0) and 16xRC bundles. The bundles show 
step-by-step structural transition. We define transition 
pressures P t start and P £ end to denote the start and end of 
the transition zone. For the 16x(23,0) Pf art = 0.5 GPa 
and P^ nd = 0.95 GPa, and for the 16xRC they are 0.45 
and 0.95 GPa, respectively. With unloading, all tubes 
remain collapsed until the return pressure P r , 0.05 GPa 
for the 16x(23,0) and 0.2 GPa the other, below which 
each tube expands to the initial structure with big hol- 
low space inside the tube. The transition pressures are al- 
most the same because they are mainly radius-dependent 
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FIG. 1. Volume vs pressure for the 16x(23,0) and 16xRC 
bundles. Vb is the volume at zero pressure. The bundles 
show step-by- step structural transition, shown in filled circles 
and squares. With unloading the pressure (open circles and 
squares), the bundles remain collapsed until the return pres- 
sure below which the bundles expand to the initial structure. 
Insects (a)-(d): Snapshots of the 16 x (23,0) bundle during 
the pressure loading up to 0.5, 0.7, 0.8, and 1.0 GPa, re- 
spectively, (e): Collapsed structure is remained when the 
16x(23,0) bundle is unloaded down to 0.05 GPa. (f): Col- 
lapsed 16xRC bundle at 0.2 GPa. (g): Larger tubes as the 
16 x (40,0) remain collapsed even under zero pressure. Here 
T = 300 K. 



[221 ], while the return pressures are different due to the 
radius inhomogeneity. For bundles with large tube ra- 
dius, for example the 16 x (40,0) whose V-P curve is not 



shown here. 
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P r < 0. Therefore it is possible to get collapsed CNT 
bundles under atmospheric pressure. The clear experi- 
mental evidence of collapsed nanotubes was very recently 
reported by the Windle group [11] from the observation 
of the "dog-bone" cross section of double- and triple- wall 
nanotubes with equivalent diameter larger than ~5 nm. 
It presents a possible way to improve the fibers during 
the direct spinning process, and it also will be of great 
importance to find a similar way for the spinning out 
of CNT forests. Insects (a)-(e) in [1] show the structural 
changes during pressure loading and the unloaded struc- 
ture for the 16 x (23,0) bundle. The unloaded 16xRC 
and 16 x (40,0) bundles with collapsed tubes are shown 
in [I] as insects (f) and (g). The herringbone structure is 
obtained due to our large compression rate [H[ . 

As a model study, we fix one of the tubes and pull an- 
other as far as possible under an external pulling force F, 
along which we define the positive direction. We find, due 
to the strong commensurability of the 16 x (23,0) bundle, 
the tube under pulling does not slide until F is larger 
than the so-called depinning force / statlc (static friction). 
[2] shows the sliding velocity at different pulling forces, 
each value extracted from a simulation longer than 100 
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FIG. 2. Sliding velocity of the sliding tube as function of 
the pulling force F under different loadings for the 16 x (23,0) 
bundle. Speeds are all abstract values, and those smaller than 
0.01 A/ps are of the same magnitude of the error by averaging 
within finite time, and thus correspond to zero velocity. Below 
the gray zone the tube stays static, and above it the tube 
slides at a speed plateau of 2600^3000 m/s, depending on 
the pressure. 



ps. Velocities smaller than 0.01 A/ps, below the gray 
horizontal zone inO are considered to be at rest because 
such values, either positive or negative, are within the 
systematic error due to the finite time average. The tube 
only slides when F goes beyond the pressure-dependent 
depinning force which will be discussed below, and the 
tube speeds up, crosses the gray zone quickly, and fi- 
nally slides at a speed plateau of 2600^3000 m/s due to 
phonon excitations. The speed plateau and its mecha- 
nism have been reported very recently [29|, [3(j , and are 
beyond the scope of current study. Obviously, such huge 
sliding speed is unreal, and as a result the bundle breaks 
due to the sliding. The sliding-induced breakage has al- 
ready been observed as an abrupt diameter change in the 
tensile- loading experiment |3|, and should be ubiquitous 
in fibers. Those speeds inside the gray zone in [2] are ac- 
tually unstable as we observed from simulation that the 
tube stops and slides intermittently. More interesting 
and important is / statlc changes with loading. Before the 
transition, it is 0.1 me v/A per atom, and is increased 
to four times when collapse happens. With unloading, 
ystatic var i es between 0.4 and 0.45 meV/A, as the case 
under 0.05 GPa shown in[2j 

With considering the change of the cross-sectional area 
of different structures one can estimate the increase of the 
tensile strength cr, 
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where the summation is over all atoms of the tube under 
pulling and A is the total cross-sectional area divided by 
the tube number in the simulation box. This equation is 
valid because within the breaking strain (less than 10%) 
of CNT bundles Q individual tubes show nearly longi- 
tudinal bond elongation rather than bond breaking and 
rotation and thus still have nearly the TPa stiffness [Ej]. 
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FIG. 3. Partial pair distribution function g(f) before nan- 
otube collapse (0.5 GPa), after the collapse (1.0 GPa), and 
after the unloading (0.05 GPa), of the 16 x (23,0) bundle. 



It means that, our assumption, the summation of pulling 
forces is identical to the external load, is correct. When 
the load becomes larger than a, it can not be totally 
transferred to neighboring tubes as the sliding happens. 
Before the transition, i.e. at 0.5 GPa, A = 378 A 2 and 
(7 = 31.2 MPa for present simulation box length about 
34 A. To reach the strength of 10 GPa, the bundle should 
be 1090 nm long, which agrees very well with the exper- 
iment jjl and of the same order as a recent theoretical 
study [31]. Now with pre-pressing to fully collapsed and 
with unloading down to 0.05 GPa, A reduces to 232 A 2 , 
a is improved by 4 x 1.63 = 6.52 due to the quadruple in- 
crease of / statlc and the ratio 1.63 between area changes. 
Thus the 10 GPa tensile strength can be achieved by a 
bundle longer than 1090/6.52 « 167 nm because all the 
load can be transferred between tubes. 

One reason for the increase of friction is the larger 
tube-tube contact area for collapsed tubes. However, it 
is not enough to make the quadruple friction increase 
because the contact goes up much less than twice. To 
investigate the structure changes, we define the partial 
pair distribution function (PPDF), 
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where a and (3 denote different tubes with atom numbers 
N a and Np, respectively, N to t the total atom number, V 
the volume of simulation box, i and j the carbon atoms 
of tube a and /?, and r$ — fj the displacement. The re- 
sult, averaged among all tubes, is shown in[3l The first 
distribution peak does not show clear changes before and 
after the structure transition, i.e., from 0.5 to 1.0 GPa, 
because it reflects the averaged tube-tube distance about 
3.5 A by using the current intertube potential. There are 
changes in the second peak, as its position shifts from r 
= 4.45 A to 4.25 A and the strength drops greatly af- 
ter the transition. Even with unloading, such changes 
still exist except that all peaks become wider and the 
total PPDF curve smoother. The third one around r = 
4.75 A at 1.0 GPa is a new peak and an evidence of AB 
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FIG. 4. Vibrational density of states of the 16 x (23,0) bundle. 
Peaks numbered as 1, 2 and 3, for example, shift to higher 
energy after the transition (1.0 GPa) and after unloading to 
0.05 GPa. These peaks are graphite-like modes, so the shift is 
larger at 1.0 GPa because the tube layers are more flattened. 
The energies are higher than the experimental values (about 
1600 cm -1 ) due to the current classical atomic potentials. 



stacking between intertube graphite-like layers. Let us 
consider graphite structure and compare to the PPDF 
at 1.0 GPa. In the AB stacking with spacing h , the 
first-, second-, third-, and fourth- nearest distances from 
atoms of one layer to those of the other, are /i, \/h 2 + a 2 , 
Vh 2 + 3a 2 , and \]h 2 + 4a 2 , respectively, where a = 1.42 
A is C-C bond length. The number ratio between them 
is 1:9:6:9. Taking h = 3.5 A from our simulation, one can 
find the second-, third- and fourth- nearest distances cor- 
respond exactly to the three PPDF peaks, respectively, 
with proper strengths. The nearest distance h is so close 
to the second one, with small number ratio as well, that 
the h -peak is overlapped as shown in [3j Here we just 
plot the PPDF for r < 7 A in order to make clear the 
structural changes. If we extend the range to r > 10 
A, each PPDF curve reaches a constant that depends on 
the atom density. For example, the PPDF constant at 
1.0 GPa is almost twice of that at 0.5 GPa and slightly 
greater than the unloaded structure at 0.05 GPa. 

The graphite-like structure of the collapsed CNT bun- 
dles can be also verified from calculations of the vibra- 
tional density of states (VDOS), 



D z (u) 



-iut 



(v z (t)v z (0)) dt, 



(3) 



where D z (ou) denotes the VDOS along the z axis and 
v z (t) the velocity of atoms along z. [H shows the graphite- 
like modes before and after the transition and after un- 
loading to 0.05 GPa. Three normal modes are labeled to 
show the energy shift, which correspond to the G band 
shift that can be observed in experiments [H, [33[ . The 
shift is larger at 1.0 GPa because the tubes are more 
graphite-like, as shown in [U After the unloading, it is 
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FIG. 5. Sliding friction under different pressures for the 
16xRC bundle. Frictions are extracted from Newton's equa- 
tion, see the text. Frictional force goes up by a factor of 1.5 
when the transition happens and also after the unloading. 



still detectable thus Raman scattering might be an effi- 
cient way to detect the tube collapse. 

However, for the worst case of commensur ability, the 
16xRC bundle, the intertube frictional force smaller 
than 0.002 me v/A per carbon atom is about two orders 
of magnitude smaller, in agreement with experiments 
[Til . Furthermore, in our pulling simulation, it is 

the sliding friction rather than the depinning static force. 
However, we still find the pre-pressing also increases the 
sliding friction between tubes. The tube under pulling 
starts to move from the optimized structure when F = 
0.005 meV/A is applied, either in or out along the tube 
axis z. The friction / is extracted from Newton's equa- 
tion when the sliding speed of that tube increases from 
zero to around 1 A/ps. For example, at 0.1 GPa, the 
speed goes up linearly from zero to 0.8 A/ps in 300 ps, 
corresponding to an acceleration of a = 0.00267 A/ps 2 . 
The sliding friction is / = ma — F = —0.00168 me v/A, 
opposite to the pulling direction. However, after unload- 
ing to 0.2 GPa, it takes more than 420 ps to reach the 
speed of 0.8 A/ps from zero, and the friction is extracted 
as -0.00263 meV/A. In[5l we plot the friction forces under 
different pressures where error bars indicate the friction 
fluctuation among several pulling simulations. Before the 
transition, frictional forces are all around 0.0017 me v/A. 
It goes up greatly to 0.0026 meV/A, by a factor of 1.5, 
when the transition happens. The friction value main- 
tains after the transition and also after the unloading. 
Taking into account the cross- sectional area changes, the 
strength of a bundle composed of tubes with random chi- 
rality is almost tripled. 



CONCLUSION 

In summary, we have investigated the mechanical ax- 
ial translation between nanotubes within CNT bundles. 
With a treatment of pre-pressing, the collapsed and com- 
pact structure remained under a pressure as small as 




FIG. 6. Tube arrangement of the 16xRC bundle in the simu- 
lation box. Primed tubes are copies of those unprimed. Tube 
1 is fixed to move while tube 3' is driven under an external 
pulling force acting on each atom. All other tubes are free to 
move. 



TABLE I. Tube chirality, radius (R), length (L), and number 
of atoms (N) in present simulation. Note that L changes 
slightly in simulation due to the barostat. 

Tube Chirality R (A) L (A) N 



16x(40,0): 


1-8 


(40,0) 


15.66 


21.30 


800 


16x(23,0): 


1-8 


(23,0) 


9.00 


34.08 


736 


16xRC: 


1 


(14,14) 


9.49 


104.00 


2240 


16xRC: 


2 


(16,12) 


9.52 


104.00 


2368 


16xRC: 


3 


(18,8) 


9.03 


104.00 


2128 


16xRC: 


4 


(20,6) 


9.23 


104.00 


2224 


16xRC: 


5 


(22,2) 


9.03 


104.00 


2128 


16xRC: 


6 


(22,4) 


9.49 


104.00 


2352 


16xRC: 


7 


(23,0) 


9.00 


104.00 


2208 


16xRC: 


8 


(24,0) 


9.39 


104.00 


2304 



comparable to atmospheric pressure. After the collapse 
of tubes, not only the cross-sectional area decreases, but 
the tube-tube frictional force also goes up, especially by 
a factor of about 4 if the tubes are commensurate. Our 
study suggests a new and efficient way to reinforce the 
strength of CNT fibers, and actually we have already 
been on the avenue of experimental studies. 



METHODS 

Periodic boundary conditions are used in all three di- 
mensions. The tubes are initially assembled in hexag- 
onal symmetry, with indexing numbers for the 16xRC 
bundle shown in [6] where primed tubes have the same 
chiralities as those unprimed. We list the chirality, ra- 
dius, length, and atom number of each tube in U The 
length is sufficiently long to effectively reflect the energy 
transfer between tubes and to avoid the size effect where 
self-diffusion can be caused for short tubes due to the 
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small energy variation between neighboring tubes. One 
constraint is used to fix the center of mass of tube 1 (see 
[6]), while tube 3' is driven to move along tube axis z un- 
der an external pulling force F acting on each atom. We 
extract the static frictional force / statlc 5 if possible, by 
assigning the critical pulling force below which no tube 
sliding happens and slides otherwise. The sliding fric- 
tional force / is extracted from the displacement-time 
curve of tube 3' by using Newton's equation F + f = ma, 
m being the mass of carbon atom and a the overall ac- 
celeration along z. Temperature T = 300 K and pres- 
sure are globally controlled by the Berendsen's algorithm 
[35^ We set along z zero pressure, P z = 0, while the 
cross-sectional pressure P x = P y varies between and 
1.2 GPa. The homogeneous isothermal compressibility 
in the cross-sectional plane is chosen to be one order of 
magnitude larger than that along z. The intratube C- 
C covalent bonds are described by the reactive empir- 
ical bond-order potential [36], while the intertube van 
der Waals interactions and those between intratube non- 
neighbor atoms by a Lennard- Jones potential (3?| , as our 
previous study has used [22]. 
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